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Abstract 
We  examine  the  rate  constants  of  chemical  exchange  reactions  involving 
relatively  light  atoms.   Such  reactions  have  specific  quantiim  characteristics, 
so  that  the  usual  classical  statistical  theory  of  chemical  reaction  rates  is 
not  applicable.   We  show  how  one  can  set  about  calculating  reaction  rates 
for  such  processes.   At  present  the  limitation  to  quantitative  calculations 
lies  in  our  inadequate  knowledge  of  interatomic  potential  energy  functions . 
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1.    Introduction 

We  shall  discuss  the  rate  constant  for  a  simple  bimolecular  reaction  - 

(1)  AB  +  C  ->  A  +  BC, 

where  A,B,  C  are  atoms  or  molecules.   The  rate  constant  for  this  reaction  vill 
in  general  be  given  by  an  expression  of  the  form 

-E   /kT 

(2)  k^jT)   =  P^^Ze  -^     , 

where  Z  is  the  number  of  collisions  of  AB  and  C  per  second  ('-'10    ,  or  less 
depending  on  density);  it  is  defined  in  terms  of  an  effective  molecular  radius 
and  thus  is  clearly  not  a  very  precise  concept.   E  f    (  ^   5-50  Kcal/mole,  i.e., 
0.2-2  ev)  is  the  'activation  energy',  smd  P„   ( '^  10   -l)  a  'steric  factor'. 
Both  P„  and  Z  will  be  mildly  temperature  dependent,  but  the  principal  -  and 
very  striking  -  temperature  variation  of  the  rate  constant  arises  from  the 
activation  energy.   The  problem  is  thus  to  explain  the  character  of  the  ex- 
pression (2)  for  the  rate  constant. 

The  standard  interpretation  of  (2)  is  that  there  must  exist  an  'activated 
state'  (ABC),  whose  energy  is  E    above  that  of  AB  +  C  at  infinite  separation. 

Q.CT« 

We  say  that  this  state  must  be  formed  if  there  is  to  be  a  reaction  (l).   Once 
the  state  has  been  formed  it  will  have  some  definite  (statistical)  life  time, 
after  which  it  can  presiomably  break  up  into  AB  +  C  or  into  A  +  BC,  both  of 
which  will  in  general  be  energetically  accessible  from  (ABC). 
Thus  the  reaction  (l)  should  really  be  written 

(3)  AB  +  C  ^->   (ABC) 

Jo 

(k)  (ABC)  r— >  AB  +  C 

oj 

(5)  (ABC)  T_->  A  +  BC 

^fj 
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where  we  use  the  index  o  for  the  state  AB  +  C,  j^  for  (ABC),  and  f  for  A  +  EC. 

Let  us  now  ask  whether  the  interpretation  (3-5)  Is  in  fact  a  necessary- 
consequence  of  the  expression  (2)  for  the  rate.   The  answer  seems  quite  de- 
finitely to  be  yes.   For,  consider  how  an  expression  like  (2)  is  arrived  at. 
Suppose  that  a  molecule  AB  and  an  atom  C  in  given  internal  states  collide  with 
relative  velocity  v^.   Then  we  calculate  the  transition  probability 
k„  (v  )  =  V  ^q(v  ),  where  6^  is  the  reaction  cross  section.   From  this 
we  get  the  temperature-averaged  transition  rate  k   (T)  by  averaging  k   (v  ) 

over  a  Maxwell  distribution  in  v  : 

o 


(6)  ^fo^^^  =(const.)  X   k^^(v^)exp 


-M  v^/2kT 
e  o 


2  , 

V  dv  . 
o   o 


Now  the  question  is,  how  can  we  get  an  expression  of  type  (2)  from  (6); 

and  the  answer  is  that  if  k_  (v  )  or  C".  (v  )  are  reasonable  functions,  then 

lo  0      fo  o  ' 

essentially  the  only  possibility  is  that  these  fiinctions  must  vanish  for 


(7)  v^  <  y  2E   ,/M      . 

o     '    act  e 

This  clearly  indicates  that  E    represents  some  kind  of  threshold  below  which 
no  reaction  can  take  place .  The  threshold  presumably  refers  to  some  quantum 
state  of  the  three  atoms  A,  B,  C  (since  it  is  not  an  eigenstate  of  the  separate 
systems),  and  this  state  is  customarily  called  the  activated  state  (ABC).   It  is 
clear  that  while  (ABC)  is  a  perfectly  definite  state  (as  this  argument  shows), 
yet  its  life  time  is  not  very  long,  because  it  is  not  observed  directly  by  normal 
chemical  techniques. 

There  are  a  number  of  steps  involved  in  calculating  the  rate  constant 
of  a  process  of  type  (l]  «iamely: 

(i)   Analyze  the  character  of  the  activated  state  (ABC), 
(ii)  Calculate  the  transition  rate  k.  (v  )  of  the  process  (3). 
(iii)  Make  some  kind  of  estimate  of  the  probability  p   of  the  activated  state 
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breaking  up  in  the  manner  (5)  rather  than  (k) , 

(8)  P,.  =lc^.  /Ufj-^oj)- 

(iv)  Finally,  carry  out  the  temperature  average  (6). 

Of  these  steps,  (iv)  is  a  trivial  technical  matter  once  the  rest  of  the 
work  has  been  done;  (iii)  is  not  difficult  in  principle,  once  we  have  solved 
step  (ii)  -  though  in  practice  we  shall  use  a  short-cut  (cf.  Section  5).   Ihis 
brings  us  to  the  major  problems  (i)  and  (ii). 

Ihe  activated  state  has  been  analyzed  by  Pelzer  and  Wigner,  London,  Polanyi, 

Pi  2l 
Eyring  and  others  i-  '  -!  .   These  authors  have  shown  that  if  one  regards  the 

nuclei  as  at  rest  in  the  field  of  the  orbital  electrons,  and  calculates  the 

total  energy  of  the  electrons  quantiim-mechanically  for  various  positions  of 

the  nuclei  (Born-Oppenheimer  adiabatic  separation),  then  one  obtains  the 

potential  energy  diagram  of  Figure  1  with  a  saddle  point  between  the  asymptotic 

N  1 
valleys  o  (=  AB  +  C)  and  f  (=  A  +  BC). '  The  state  j    =  (ABC)  is  regarded  classi- 
cally as  a  state  located  on  or  about  the  saddle  point. 

The  problem  was  taken  up  quantum-mechanically  by  Bauer  and  Wu  ^  J ,  who 
followed  a  suggestion  of  Mott  and  Massey^ J .   The  state  (ABC)  is  regarded  as 
vibrationally  excited,  which  is  reasonable  in  view  of  its  character  involving 
nuclear  motion,  and  of  the  order  of  magnitude  of  the  energy  of  excitation, 

N.l 

In  Figure  1,  as  in  most  of  the  classical  calculations  and  in  the  present  con- 
siderations, we  have  assumed  that  the  three  atoms  A,  B,  C  lie  in  one  straight 
line  throughout  their  collision.   If  one  tried  to  consider  glancing  collisions 
and  rotational  motion,  the  'potential  energy  surface '  of  Figure  1  and  Section 
3  would  have  to  be  replaced  by  a  six-dimensional  hypersurf ace .   Some  effects 
that  would  occur  if  one  tried  to  discuss  rotational  motion  are  mentioned  in 
Section  8. 
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AB 


FIGURE  I 


Potential  energy  surface  of  a  collinear  system  A-B-C,  in  terms  of  the  coordin- 
ates r   and  r .  The  arrows  indicate  the  approach  of  C  towards  AB  and  the 

AB      BC 

final  separation  of  A  from  BC.  B*  is  the  saddle  point  of  the  surface. 


N.2 
E   .  .   Now,  evidently  a  state  of  this  kind  has  to  be  treated  quantum- 

SlCC 

mechanically  provided  its  vibrational  quantum  number  is  low  (of  order  1, 
not  10) .   In  fact  for  interactions  of  light  atoms  the  vibrational  quantum 
number  is  low^-^-',  and  thus  we  have  an  immediate  reason  why  the  statistical 

Ti  2I 

treatments  ^  '  -^ ,  which  treat  the  activated  state  classically,  cannot  be  re- 
garded as  a  satisfactory  approximation.   In  fact  there  are  very  serious 
difficulties  in  attempting  a  detailed  quantum-mechanical  calculation,  but 
now  one  can  at  least  analyze  them.  Let  me  remark  that  at  present  the  essential 
limitation  lies  in  the  inadequacy  of  presently  available  quantum-mechanical 
calculations  of  potential  energy  surfaces.  Tb.is   will  be  discussed  in  more 
detail  in  Section  3. 

The   final  point  is  (ii),  the  calculation  of  the  transition  rate  ^^^('^q)' 
Here  there  is  very  definite  progress  to  report  over  the  early  work  ^  J .   If  we 
treat  the  activated  state  as  vibrationally  excited,  then  recent  calculations'--^'  '  -■ 
enable  us  in  principle  to  make  proper  quantum-mechnical  calculations  of  the  rate 
constant  which  are  essentially  limited  only  by  our  inadequate  knowledge  of  the 
potential  energy  surfaces. 

It  is  clear  that  there  is  a  strong  qualitative  similarity  between  the 
present  problem  and  Bohr's  compound  nucleus  theory,  which  describes  the  inter- 
action of  nucleons  with  relatively  heavy  nuclei.  However,  there  is  an  essential 
quantitative  difference  in  the  quantiim  character  of  the  activated  state  for  groups 
of  light  atoms,  and  of  the  compound  nucleus  containing  a  fairly  large  niunber  of 
nucleons  (N  >  25).   The  compound  nucleus  has  very  many  closely  spaced  excited 

N.2 

A  characteristic  energy  for  rotational  degrees  of  freedom  is  '^.05  ev, 

-^  0.2  -  0.5  ev  for  vibration;   '^  5  ev  for  electronic  excitation.  Energies 
of  dissociation  (a  vibrational  type  process)  are  of  order  5  ev,  while  ioniz- 
ation potentials  are   ~  10  ev. 
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states,  any  or  all  of  which  can  be  occupied,  so  that  the  ensemble  can  be 
treated  statistically  or  'classically'.   In  the  present  case  the  activated 

state  has  very  few  quantum  states,  and  must  therefore  be  treated  by  quantum 

N.3 
mechanics  . 


2.   Formulation 

To  discuss  the  problem  further,  we  must  give  a  proper  formulation.   In 
the  statistical  theory  it  is  customary  to  use  the  coordinates  r^  and  Ic"  p,  but 
in  a  quantum-mechanical  calculation  we  wish  to  eliminate  the  motion  of  the 
center  of  mass,  and  also  to  write  the  kinetic  energy  as  the  sum  of  squares  of 
momenta,  without  any  cross  terms.  To  achieve  this,  we  use  the  coordinates^  -' 

^  =   (^A  ^A  ^  ^  ^B  ^  "C  ^C^/(^A  ^  ^  ^  "C^ 


(9) 


for  which  we  write  the  quantum  hamiltonian  as 


(10) 


(a) 

i      (b) 

(c) 


H  =   H%  H°+  H°  +  V'(p,?) 

H°  =   (-b^/2M.)^  +  V°(^);  M.  =  M^Mg/(M^+Mg) 


H°  =   (-ti^/2M  )^  +  V°(r);  M 
e       '   e  r    e      e 


^c(v^)/(w\) 


N.5 

If  one  considers  glaincing  collisions,  i.e.,  rotational  fine  structure,  there 

will  be  many  more  states  than  in  the  present  discussion.  However,  the  selection 

rules  for  rotational  transitions  (cf .  Section  8)  indicate  that  for  a  single 

tretnsition  process  the  situation  will  be  essentially  as  outlined  here,  although 

the  potential  energy  curves  will  be  slightly  different  for  different  values  of 

angular  momentvun. 
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H°  =  H^(R)  gives  only  the  motion  of  the  center  of  mass,  which  does  not  affect 
the  cross  section.   It  is  omitted,  and  thus  conservation  of  linear  momentum 
is  automatically  ensured  and  does  not  have  to  be  considered  in  the  following. 
Of  the  two  remaining  coordinates,  'p   describes  the  internal  motion  of  the  mole- 
cule AB,  ('i'  for  internal),  while  r  represents  the  external  ('e')  motion  of 
C  relative  to  AB. 

Now  we  wish  to  solve  the  quantum- mechanical  problem 

(11)  (H°  +  H°  +  V  -  E)  fjt,r)   =0;    n  -  o^,i,f. 


In  fact,  we  shall  solve  the  problem  by  perturbation  theory,  basing  the  dis- 
cussion on  eigenfxinctions  of  the  uncoupled  hamiltonian  (H.  +  H  ) '^-' , 


/ 


(12)  ^ 


(a)  (H°  +  H°  -  E)  $  jt,T)    =  0. 

(b)  $j^(p,?)   =  0^(p)iir^(r) 

(c)  (H°  -  e^)0^(p)  =  0 


(d)  r(r)  =  "^    XA\r)?  (cos  Q) /t 


(e) 


,2/,  2 


2.„o, 


d^/dr""  +i(/+l)/r^-(2M  /fe'')V°(r)+  k 


X,(k  r)  =  0 


N.4 


V 


(f)  JL^ik^r)  ^  cos(k^r  +  &„  ^  )  as  r  -> 


Conservation  of  energy  and  angular  momentum  between  two  states  n,  n'  gives 
the  conditions 


(13) 


(1^) 


E  =  (ti^  k^  /  2M  )  +  e    =  (fe^  k^,/2M  )  +  €  ,  , 
^    n  '   e      n     ^    n"   e'     n'  ' 


A,=  ^n' 


N-^  Note  that  in  writing  (l2e)  we  have  assumed  V°(r)  =  V°( |r | )  =  V°(r).   If 

this  is  not  so,  one  should  write  V  ("r)  =  T"  V   (r)P  (cosO  ),  and  proceed 
'  e  _     e,m^  '  ra^    n  '  ^ 


m=o 


from  there. 
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Ihe  transition  probability  per  unit  time  is  given  by 


(15)     K<  ^'^ 


n'.n     S~ 


(  ^n"  ^'  ^n^l   '^P/'^^(n') 


where  dp/dE/  ,.  gives  the  density  of  final  states  per  unit  energy  rsinge. 

We  are  here  setting  up  the  problem  in  terms  of  quaintum-mechanical  pertur- 
bation theory.   The  use  of  quantum  theory  here  is  of  course  distinct  from  its 
use  in  describing  the  activated  state,  but  in  fact  a  quantiom  treatment  is  again 
necessary'-  -■ ;  this  point  is  discussed  in  detail  in  Section  6.  One  could  of 
course  solve  the  classicail  problem  corresponding  to  (9)  -  (l5)j  although  a 
considerable  reformulation  would  be  necessary. 

Regarding  the  use  of  perturbation  theory,  this  is  essentially  unavoidable. 
Even  in  classical  mechanics,  the  three-body  problem  cannot  be  solved  exactly, 
in  general,  and  thus  we  can  hardly  expect  to  do  any  better  in  quantum  theory. 
Some  earlier  quantum  calculations  of  the  excitation  of  molecular  vibration 
by  collisions  talk  about  using  'distorted  wave'  or  'perturbed  stationary 
state'  methods;  these  are  not  exact,  but  somewhat  better  in  principle  than 
straight  perturbation  theoiy.  However,  when  it  comes  to  the  point,  everybody 
uses  some  form  of  perturbation  theory.  Ihere  is  one  serious  difficulty  in 
principle  which  can  be  important  in  certain  energy  ranges.  The  scattering 
matrix  is  not  unitary  in  perturbation  theory,  so  that  the  calculated  trans- 
ition probabilities  do  not  satisfy  flux  conservation  unless  K^j^.V  ^^)  | 

N.5 
is  small  in  some  appropriate  sense. 

This  problem  is  not  likely  to  arise  in  the  present  application  at  reason- 
able energies,  but  it  will  enter  into  discussions  of  three-body  recombination 
where  we  are  dealing  with  relatively  small  energy  transfer  and  consequently 
large  transition  probabilities.  Work  on  this  problem  is  being  done  (in  different 
contexts)  by  J,  Keck  at  AVCO  and  by  M.  Salkoff  at  New  York  University. 
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Before  we  a,o   any  further,  we  have  to  examine  the  potential  energy  in 
detail.   It  is  known  ^  -'  that  the  cross  section  is  an  exceedingly  sensitive 
fimction  of  the  potential,  and  thus  we  shall  examine  the  problem  of  obtaining 
an  expression  for  the  potential  energy,  and  of  separating  V',from  V.  and  V  . 
This  problem  is  discussed  in  Section  ^.      In  the  following  sections  we  shall 
examine  the  kind  of  model  for  the  activated  state  *  used  in  the  statistical 
theory  L  '  J  and  in  the  earlier  quantum  mechanical  work  '-^-' . 

3.   DifficTilties  in  calculating  the  potential  energy 

It  is  quite  clear  in  principle  how  one  sets  about  calculating  the 
potential  energy  function  which  determines  the  motion  of  the  nuclei  in  the 

' smeared-out '  field  of  the  orbital  electrons.  One  makes  a  quantxjm-mechanical 
calculation  of  the  energy  of  the  electrons  moving  in  the  field  of  the  sta- 
tionary nuclei.  By  making  this  calcvilation  for  a  variety  of  (fixed)  positions 
of  the  nuclei,  one  obtains  an  energy  function  which  determines  the  interaction 
of  the  heavy  particles  with  one  another.   In  particular,  if  one  considers  the 
three  nuclei  A,  B,  C  as  moving  in  a  straight  line,  then  the  potential  energy 
ftinction  will  be  the  usual  'potential  energy  surface'   '  '    . 

The   difficulties  in  carrying  out  1±Lis  program  are  'merely'  technical,  but 
they  are  exceedingly  formidable  on  account  of  the  large  niimber  of  electrons 
involved.  For  instance,  for  the  cases 

If  we  consider  several  electronic  states  of  the  system  as  a  whole,  then 
each  will  have  its  own  potential  energy  function  or  surface.   We  just  consider 
the  lowest  of  this  manifold  of  surf aces, disregarding  electrpnic  excitation. 
Ihis  is  justified  by  the  adiabatic  hypothesis;  in  practice  it  is  very  satis- 
factory if  we  restrict  ourselves  to  temperatures  below  about  5OOO  K. 


(16) 


^ 


(a)   H   +  F  — >  HF  +   H 


(b)  N   +  0  — >  NO  +  N 


there  are  respectively  11  and  22  orbital  electrons,  giving  three  times  this 
number  of  degrees  of  freedom;  even  with  the  simplification  of  considering 
only  electrons  in  the  outermost  unfilled  sub- shell  in  each  case,  we  still 
have  9  and  11  electrons  respectively.   An  actual  solution  of  the  wave  equation 
for  such  a  system  by  either  analytic  or  niimerical  means  seems  to  be  out  of  the 
question  with  presently  available  techniques,  and  thus  one  is  reduced  to  using 

•  +•   T  +  V  •    [lO^lll 
variational  techniques  •-    -*  . 

The  accuracy  required  is  of  the  order  of  thermal  energies,  i.e., -03  ev, 
as  compared  with  molecular  binding  energies  of  the  order  of  3  ev.  Such 
accuracy  is  exceedingly  difficult  to  obtain;  for  the  classic  case  of  the 
ground  state  of  the  hydrogen  molecule  ^  -" ,  which  involves  only  two  electrons, 
such  orders  of  accuracy  were  only  obtained  with  wave  functions  involving  about 
10  variational  parameters.  This  may  well  be  a  somewhat  unfavorable  example, 
and  more  recently  developed  techniques  may  cut  down  the  number  of  variational 
parameters  required  per  electron,  but  the  problem  is  still  very  difficult. 
I  believe  that  with  modern  high-speed  computing  equipment  as  fast  as  the  IBM 
704  or  Remington  Rand  llOJ,  it  woiild  be  possible  to  attempt  the  calcvilation  of 
potential  energy  functions  to  this  order  of  accxiracy. 

Thus  the  position  is  clear.  At  present  it  is  not  feasible  to  attempt  to 
make  quantitative  calciilations  of  reaction  rates  for  even  the  simplest  processes 
because  we  lack  the  queintitative  knowledge  of  the  potential  energy  functions 
necessary  to  ensure  meaningful  answers.  However,  we  can  give  a  qualitative 
analysis  of  the  problem,  discussing  significant  features  and  needed  improve- 
ments, with  the  aim  of  showing  how  one  can  set  about  doing  these  calculations 
in  principle. 
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k.       The   activated  state  In  the  classical  statistical  theory'-/'  -* 

This  is  a  brief  review  of  the  usual  theory.   In  Figure  1  we  have 

N.l 
the  potential  energy  surface     in  its  customary  form,  in  terms  of  the 

coordinates  r._  and  r_,^.   We  consider  the  system  moving  along  the  initial 

Ad       dO 

* 
track  0-,  Op  up  to  the  saddle  point  B   where  it  stays  for  a  relatively  long  tune 

because  most  of  its  energy  has  been  converted  from  kinetic  to  potential  energy, 

so  that  its  velocity  is  very  low,  and  hence  it  takes  a  long  time  to  travel 

aroxind  on  the  saddle  point.  The  region  of  the  saddle  point  represents  the 

activated  state  ^,   and  its  height  above  the  valley  floor  at  infinity  (O, )  is 

the  classiceuL  activation  energy.  The  system  stays  in  the  activated  state 

sufficiently  long  for  it  to  'forget'  how  it  got  there;  eventually  it  goes 

either  into  the  final  state  f  or  into  the  initial  state  o. 

To  calciilate  the  rate,  one  defines  a  mean  collision  frequency  Z  =  Z(T) 

which  indicates  the  number  of  times  an  AB  molecule  and  a  C  atom  come  s\iff iciently 

close  to  one  another  (say  to  Op)  for  there  to  be  a  reaction.  Once  it  gets  into 

this  configuration, then  the  probability  of  the  system  getting  up  to  the  saddle 

will  be  proportional  to  the  relative  volume  of  phase  space,  or  to  the  partition 

function  of  the  system  'in'  the  saddle  relative  to  that  of  the  system  along  the 

O-asymptote.   If  the  partition  function  for  the  system  in  state  n  (=  o,j,f)  is 

f^(T),  then 

(IT)         vf^^^^Cr)     =    fj(T)/f^(T), 

while 

(18)  vf^^^'W     =     P^f^^(T)  vll^''iT), 

It  is  customary  to  write  p-.   (T)  as  simply  an  adjustable  K  ~  p  ;  this  was 
also  done  in  the  earlier  quantum  calculation'--'. 

Iliis  is  a  brief  outline  of  the  theory.  Many  applications  of  it  have  been 
made;  sometimes  (2)  is  replaced  by  the  equiveLLent  expression 
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f 


(a)  k^^(T)  =  Z(T)  exp 


act 


A^ 


(19)  <     (b)  F^^^  =  E^^^  -  T  S^^^ 


(c)  exp(S^^^/k)  =  P^^. 

F  .  and  S  .  are  called  free  energy  and  entropy  of  activation.  This  term- 
inology  is  used  in  the  'absolute  reaction  rate  theory',  a  version  of  the 
theory  which  describes  reaction  rates  in  terms  of  equilibritmi  thermodynamic 
concepts • 

5.   Quantum-mechaMcal  model  of  the  activated  state*-  ^ 

First  of  all,  let  us  redraw  the  potential  energy  surface  of  Figure  1  in 
terms  of  the  coordinates  (p,r)  of  (9).  Hiis  is  shown  in  Figure  2.  To  attempt 
to  solve  the  motion  of  the  system  along  a  track  0,  Op  B  Fg  F^  would  be  very 
difficult  even  if  we  knew  the  orbit  accurately,  which  we  don't.  Thus  we  make 
an  important  simplifying  asstjmption  to  enable  us  to  carry  out  the  calculation. 
Assumption: 

We  replace  the  actual  system  by  a  model  in  which  the  passage  of  the  system 
(AB  +  C)  during  the  approach  of  C  from  infinite  separation  to  a  configuration 
in  which  AB  and  C  are  very  close  together,  takes  place  adiabatically,  i.e., 

ri2i 

without  quantum  Jumps;  the  quantum  transition  0  ->  j  then  takes  place  suddenly '-  -■ 
in  the  immediate  neighborhood  of  the  configuration  of  close  approach. 

This  assumption  can  be  partly  justified  by  the  short  range  of  the  interaction 
potential  responsible  for  the  transition  (l). 

We  see  thus  that  the  process  of  attaining  the  activated  state  ^  falls  into 
two  parts.  First  of  all,  during  the  adiabatic  passage  of  the  system  from  in- 
finite to  small  separation,  the  zero-point  vibrational  energy  of  the  molecule  AB 

will  change  by  a  (presumably  relatively  small)  amount  E   :  this  energy  is  supplied 

act 


11a 


/>t=^«) 


AB 


FIGURE  2 


Potential  energy  surface  in  Fig.  1  in  terms  of  the  coordinates  p,   r  defined  by 
Eqs.  (9)'  The  arrows  have  the  same  meaning  as  in  Fig.  1.  The  kinetic  energy 
of  relative  motion  (in  the  coordinate  r)  is  assumed  to  be  transformed  into 
the  vibrational  motion  of  the  system  ABC  ( in  the  coordinate  p ) ,  forming  the 
excited  state  of  the  configuration  near  the  saddle  point  B*. 
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from  the  relative  kinetic  energy  of  the  system.  Tben,   when  the  atcsn  and 

molecule  are  very  close  together,  the  quantum  jump  from  the  ground  vibrational 

1 1 
state  o  to  the  excited  state  j  will  take  place  with  a  transfer  E^^^  of  energy 

frcMi  kinetic  into  potential  (vibrational)  energy.  Thus  the  total  activation 

energy  E  .  is  the  sum  of  these  two  contributions,  of  which  the  second  will  in 
ac  u 

N.7 
general  be  very  much  larger  than  the  first 

(20) 


t      It 

E^  =  E.+E. 
act      act    act 


I  •» 

E  ^  =  van  der  Waals'  energy,  ~'  .02  ev;  E    =  vibrational  energy,  ^   .2  ev. 
act  ac"C 

In  the  light  of  this  assumption,  we  can  now  break  up  the  potential  energy 
into  V°  +  V  in  such  a  way  that  a  calculation  can  be  made.  We  do  this  by  ex- 
panding V(p,r)  about  the  value  V(pj,rj)  at  some  point  (pj,  rj)  which  is  character- 
istic of  the  configuration  at  the  distance  of  closest  approach  of  AB  and  C,  and  we 
then  treat  (p-Pr)  and  (r-r,.)  as  small  L^J  *   .  Thus  the  total  potential  energy  is 

'  (a)  V(p,r)  =  V°(p,rj)  +  V°(pj,r)  +  v'(p,r) 

(21)   -J       ,  ^V  (p^,r)         Sv°(p,r  ) 

(b)  V(p,r)=  (P-Pj)   %'         -   (r-rj)   L  -—• 

k  Jo 

We  see  that  as  p-  and  r,  are  constants,  we  have  now  f orma.1  ly  carried  out 

a  separation  of  the  hamiltonian  into  additive  terms  in  p  and  in  r,  which  together 

I 

form  the  unperturbed  hamiltonian,  and  a  perturbation . V  which  couples  the  p  -  and 

In  our  earlier  calculations  L^-'  we  put  E  +  =  0;  the  existence  in  principle 
of  the  term  E  .  was  pointed  out  to  me  by  Dr.  B.  Widom. 


N  8 
*  Because  the  distances  p^  and  r^  are  comparable,  we  cannot  treat  p  and  r  in 


the  unsynmetric  way  that  is  used  in  some  problems 


[5,6,7] 
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r-motions,  This  separation  is  meaningful  in  the  sense  of  the  assumption  made 

above . 

Let  us  now  examine  the  potentials  V°(p)  and  V°(r)  more  closely.  The 
potential  V^(r)  is  obtained  in  principle  by  taking  a  section  through  the 
potential  energy  surface  of  Figure  2  keeping  p  =  Pj  =  constant,  i.e.,  roughly 
through  the  point  B  .   It  has  the  character  of  a  van  der  Waals  potential,  being 

strongly  repulsive  for  small  r  and  weakly  attractive  for  large  r.  The  potential 

„0/  \  .      ,  * 

v^tpj  IS  much  more  interesting.   If  we  take  a  section  through  the  point  B  of 

Figure  2,  keeping  r  =  Tj  =  constant,  we  get  a  double  minimum  potential  of  the 

N  Q 
sort  shown  in  Figure  3   '^. 

The  case  where  the  potential  energy  function  of  Figure  5  is  symmetrical 
with  respect  to  inversion  about  the  center  point  B   has  been  examined  rather 
thoroughly  before;  here  we  wish  to  examine  the  assymetric  case  qualitatively. 

First  of  all,  it  is  clear  that  Figure  3  is  just  one  of  a  sequence  of 
similar  sections,  for  various  values  of  r.  A  typical  one  for  a  very  much 
greater  value  of  r  than  rj  is  shown  in  Figure  4.   If  we  start  with  a  molecule 
AB  in  its  lowest  vibrational  state,  then  clearly  we  are  in  the  vibrational 
state  of  Figures  3,4.  We  now  ask  how  the  system  can  go  over  from  AB  +  C  into 
a  state  which  corresponds  to  (A  +  BC) .  Clearly  the  state  of  lowest  energy  that 
has  a  finite  probability  of  finishing  up  as  (A  +  BC)  asymptotically  is  2,  and 
the  probability  of  this  ending  as  (A  +  BC)  is 

(22)     Kc2  =  lim  ^  [  \0Jp)\^dp  /    [   \0Jp)fdp 

Now,  we  can  go  from  state  1  to  state  2  in  several  ways,  all  of  which  call 
for  energy  to  be  supplied  from  an  outside  source,  namely  the  translational 
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For  problems  of  this  general  character  the  minimum  will  always  be  double. 
For  values  of  r  less  than  Tj  the  curves  will  shrink  to  give  just  a  single  minimum, 
but  the  energy  of  any  eigenstates  of  this  sort  of  configuration  will  be  higher 
than  any  of  the  states  we  are  considering. 


15a  - 


^  P 


FIGURE  3 


The  potential  V.Cp^r^.),  which  is  a  section  of  the  potential  energy  surface  of 
Fig .  2  at  constant  r  through  the  saddle  point  B* • 
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FIGURE     4 


A  potential  V. (p,r'),   which  is  a  section  of  the  potential  energy  sxirface  of 
Fig.   2  at  constant  and  veiy  large  r  =  r'. 


Ik 


(thermal)  motion   '   .  The   possible  ways  of  getting  from  state  1  to  2  are: 
(i)   Excitation  by  an  amount  of  energy  of  (Ep-E,),  and  then  tunnel  effect  \intil 
the  wave  ftonction  has  leaked  through  the  central  barrier  to  such  an  extent  that 
the  system  is  in  an  approximation  to  the  eigenstate  2,  if  this  can  be  accomplished 
during  one  collision.   In  fact  we  know  that  this  is  not  an  important  process, 
theoretically  because  leakage  times  are  exceedingly  long  compared  with  collision 
times  and  experimentally  because  such  a  process  would  not  be  strongly  temperature 
dependent  for  an  exothermic  reaction,  and  would  give  E    =  (E^-E, )  for  an  endo- 
thermic  reaction.   This  is  not  observed  in  practice  for  most  processes, 
(ii)  Excitation  of  state  3-  We  see  from  Figure  h   that  this  state  is  not  stable 
as  r  increases,  since  it  will  presumably  degenerate  into  either  state  1  or  state 
2,  depending  on  'which  side  of  the  barrier  the  atom  B  is  on  during  the  separ- 
ation'. The  excess  energy  will  go  back  into  translation.  Thus,  once  the  state 
3  is  formed,  there  is  a  probability   Kpp  that  we  finish  with  the  system  (A  +  BC). 
(iii)  Excitation  of  some  higher  state  k,^,..   which  degenerates  into  state  2  just 
the  way  5  does.   In  fact,  at  reasonable  temperatures  (E. -E,)  »  kT,  so  that  we 
may  disregard  these  processes,  since  it  is  very  unlikely  that  the  system  will 
have  sufficient  energy. 

We  see  thus  that  an  essential  difference  between  the  present  quantum-mechanical 
model  and  the  usual  classical  one  lies  in  the  interpretation  of  the  activated 
state,  and  in  particular  of  the  activation  energy.  Classically  this  is  the 
height  of  the  saddle  point  above  the  potential  at  infinite  separation,  but  in 
the  present  quantum  theoretical  model  it  is  the  difference  between  two  appropriate 
vibrational  energy  levels. 


N.IO 

We  are  considering  an  endothermic  process.  We  could  equally  well  consider 

an  exothermic  process  by  just  inverting  the  argument. 
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Having  now  defined  our  model,  let  us  ask  whether  a  quantum  treatment  of 
this  aspect  of  the  problem  is  necessary  in  principle.  The  answer  is  quite 
definitely  yes,  because  the  wall  is  rather  small,  so  that  we  have  relatively 
few  vibrational  states  1,2,5,  ...  L-^-i  .  Quantitatively,  if  E  ,    is  the 
height  of  the  central  maximum  above  the  bottom  of  the  one-well  of  Figure 
3,  then  the  expression 

f 
(23)  (lA  fe) 


/2Mil(Eelass-^?^P)l  ^ 
is  of  order  1-2,  not  10-100,  so  that  a  quantum  treatment  is  necessary. 

6.   Quantum  characteristics  of  the  rate  calculation 

We  have  discussed  the  character  of  the  activated  state  on  the  present 
model,  and  have  seen  that  it  is  essentially  quantum- mechanical,  a  vibrationally 
excited  state.  Let  us  sketch  the  characteristic  features  of  the  process  of 
formation  of  this  state,  which  follows  closely  from  earlier  work  on  the  ex- 
citation  of  molecular  vibration  by  collision'-   -^: 

(i)  On  account  of  the  large  mass  of  the  interacting  molecules,  the  de  Broglie 
wavelength  for  the  relative  motion  is  small  compared  to  the  range  d  of  the 
intermolecular  potential  for  the  states  o  and  f : 

(24)  (t  /  M^  v^)   =  1/k^  =  ?c^  «  d;      n  =  o,f . 

Ihis  means  that  for  these  states  we  can  talk  of  a  quasi-classical  orbit  of  the 
particles,  and  presiimably  do  not  have  to  use  better  than  a  WKB  approximation 
for  the  wave  functions;  but  equally  we  have  to  consider  large  values  of  angular 
momentum,  at  least  up  to 


(25)  i-^  :&  =  k  d  fe  --  60  t. 
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This  point  was  not  properly  understood  in  the  early  workLi,  and  reflects  the 
significant  improvement  of  the  more  recent  calculations  L '  '  J . 

Note  that  condition  (2^)  does  not  apply  to  the  state  j,  because  essentially 
zero  kinetic  energy  implies  large  de  Broglie  wavelength,  regardless  of  the  mass. 
This  means  that  t,(r)  cannot  be  treated  in  the  WKB  approximation,  but,  because 
condition  (25)  is  satisfied  for  o-  and  f- states,  we  shall  still  have  to  consider 
large  values  of  angular  momenttmi  /fe. 

(ii)  At  the  low  energies  that  are  of  real  interest  here,  the  probability  of 
energy  trajisfer  is  small  because  there  are  many  vibrations  (of  frequency  ai/2n) 
during  one  collision,  whose  average  duration  is  of  order  d/v, 

(26)  CO  »  v/d. 

(lii)  To  get  any  further,  we  have  to  make  a  more  detailed  investigation  of  the 
mechanism  of  energy  transfer.   Let  us  use  our  perturbation  theory  model,  in 
which  the  coupling  is  induced  by  a  potential  V  .  Classically  we  may  speak  of 
a  definite  time -dependent  orbit  of  the  incident  particle,  giving  V  as  a 
function  of  time.  Then,  classically,  the  average  momentum  transferred  in 

ri5i 

a  collision  is  >-  -" 

f  _ 

(27)  6p  =    V  dt  21  V'   d/v 

J 

where  V^  is  a  suitable  average  of  V'(t).  and  d/v  is  the  collision  time.  Now, 
we  have 

(28)  6p  d  «  fe, 

N.ll 
which  is  an  explicit  quantxim  condition 


If  we  were  to  apply  this  argument  to  Coxilomb  interactions  (disregarding  the 
infinite  range  of  the  Coulomb  potential),  and  put  V  =  Z  e  /d,  then  condition 


(28)  would  just  become  Z  e  /ti  v  «  1, 
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Let  us  now  look  at  conditions  (i)-(iii)  overall.   They  are  clearly  all 
independent  conditions;  (i)  and  (iii)  are  quantiJin  conditions  which  involve 
Planck's  constant  explicitly;  condition  (ii)  is  the  'adiabatic  condition' 
which  holds  both  classically  and  quantally,  and  hence  does  not  contain  ti 
explicitly.  Thus,  the  characteristics  of  the  present  problem  as  displayed 
by  these  three  conditions  are  the  following: 

Condition  (i)  shows  that  large  values  of  angular  momentum  Z  b  have  to  be 
taken  into  accoixnt. 

Condition  (ii)  shows  that  the  transition  probability  will  be  small,  so  that 
a  careful  treatment  is  necessary. 

Condition  (iii)  supplements  condition  (i)  and  shows  that  a  classical  treatment 
of  the  collision  is  not  adequate. 

7-   Sketch  of  possible  methods  of  solution 

Given  some  reasonable  approximation  to  the  potential  energy  surface,  one 
can  carry  out  the  present  sort  of  rate  calculation  with  a  computing  machine 
as  fast  as  a  UNIVAC  or  IBM-65O.  The  reason  for  this  is  that  as  we  need  maybe 
60  values  of  angular  momentum   -£  fe,  it  is  necessary  to  solve  the  radial  wave 
equation  (l2e)  for  at  least  5  values  of  jL    to  have  something  to  interpolate 
between.  Ihis  has  to  be  done  for  both  initial  and  final  wave  niunbers,  so 
that  at  least  10  translational  wave  functions  have  to  be  found  for  every 
point.  Numerical  integration  of  the  wave  equation  takes  about  3  minutes  for 
one  value  of  k  and  jL   ;  writing  down  the  WKB  wave  function  will  take  about  one 
minute  per  point,  as  will  the  evaluation  of  the  treinsition  matrix  element,  once 

r7i 

the  wave  functions  are  given ^ J.  I^us  the  time  for  computing  translational 
matrix  elements  for  one  value  of  initial  energy  will  be  of  the  order  of  3O-6O 
minutes. 
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Actually  the  evaluation  of  translational  wave  functions  and  matrix  elements 
is  inherently  rather  well- suited  to  analog  computing  equipment.  Some  exploratory 
calculations  carried  out  by  Mr.  J.  P.  Krawiec  at  AVCO  show  obvious  advantages  of 
analog  over  digital  techniques.   With  what  I  understand  to  be  excellent  equipment 
and  veiy  good  programming,  the  time  required  per  matrix  element  (which  includes 
finding  the  wave  functions)  will  be  of  the  order  of  10-50  minutes. 

On  digital  equipment  one  starts  the  nimerical  integration  at  the  turning 
point  of  the  classical  motion  with  the  values  of  function  and  derivative  given 
by  the  Langer-WKB  solution;  the  exponential  decrease  of  the  wave  fxinction  for 
smaller  distances  is  approximated  by  a  WKB  wave  function.   On  analog  equipment 
one  can  do  much  better  than  this  in  principle.  The  equipment  is  set  up  to  pro- 
duce initial  and  final  wave  fvinctions  simultaneously  and  to  evaluate  the  matrix 
element  in  the  same  operation  .  The  method  of  finding  the  solution  is  to  vary 
the  ratio  of  wave  function  to  derivative  at  one  point  until  one  gets  a  solution 
that  comes  veiy  close  to  zero  as  one  integrates  backwards,  towards  the  origin. 
This  is  a  very  sensitive  test;  the  equation  (l2e)  has  two  linearly  independent 
solutions,  of  which  the  regular  one  (which  we  want)  goes  to  zero  at  the  origin, 
while  the  irregular  one  diverges.  Itus,  if  one  integrates  backwards,  even  the 
smallest  admixture  of  the  irregular  solution  will  give  a  wave  function  which 
diverges  at  the  origin.  The  sensitivity  of  the  UNIVAC  is  so  high  that  inte- 
grating backwards  always  gave  solutions  that  diverged  very  sharply  to  +  oo  ; 
on  the  analog  machine  it  has  proved  feasible  in  terms  of  time  to  search  for  a 
solution  which  comes  quite  close  to  zero  before  it  diverges. 

An  obvious  objection  to  analog  equipment  lies  in  the  accuracy  that  can 
be  obtained;  the  results  obtained  are  reproducible  to  better  than  1  /o  for 
the  matrix  element,  and  the  absolute  accuracy  of  the  matrix  elements  is  esti- 
mated as  within  5  /o.   I  feel  that  5-10  /o  error  is  about  the  maximum  that  can 
be  tolerated  in  the  matrix  element,  not  so  much  in  terms  of  the  accuracy  of  the 
final  numerical  results,  as  that  if  the  total  cross  section  is  iincertain  by  more 
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than  10-20  /o  it  will  be  very  hard  to  draw  any  conclusion  from  one's  results. 

In  addition  to  the  translational  matrix  elements,  one  also  has  to  find 
the  internal  wave  functions  0  (p)  and  their  matrix  elements.  This  is  not  a 
major  problem  because  it  just  has  to  be  done  once  for  the  whole  calculation. 

8.   The  effect  of  rotation'^ ''''^ 

In  the  present  analysis  we  have  considered  only  collinear  collisions  of 
the  three  particles  A,  B,  C.  This  is  a  simplification  which  is  almost  always 
made  in  the  statistical  calculations  L '  J   •  ,   j^  practice,  we  shall  find 
that  most  actual  collisions  occur  at  an  angle  rather  than  head-on,  and  so 
the  effect  of  rotational  motion  of  the  AB-molecizle  should  be  considered  in  any 
quantitative  calculation.  This  will  lead  to  the  following  differences: 
(i)   Conservation  of  angular  momentum  will  no  longer  be  just 

W  /„  =  I,,  , 

as  has  been  considered  in  the  preceding,  but 

(29)  L  +  J   =  L  ,  +  J  ,  , 

n    n      n'    n'  ' 

where  L  is  the  orbital  angiilar  momentum  in  the  r-raotion,   having  the  eigenvalue 

t  >/TJ~t+l) ,   while  J  is  the  angular  momentum  of  rotation  of  the  AB-molecule, 


i.e.,  in  the  p-motion.   Its  eigenvalues  will  be  fe  yj(j+l)  . 
(ii)  There  will  be  a  potential  energy  surface  for  each  possible  combination 
(JQ-Jf.)>  i.e.,  for  each  initial  and  final  state  of  rotation  in  the  p-coordinate. 
Thus  there  will  be  a  corresponding  change  in  the  condition  for  conservation  of 
energy,  (15),  in  that   e^  =  ^n^^n"^n^'  where  v  and  J  are  respectively  vibra- 
tional and  rotational  quantum  numbers. 


N.12 

Ihis  section  was  written  as  a  result  of  some  remarks  by  Dr.  S.H.  Bauer. 
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(iii)  Precisely  which  transitions  J  -j„  are  possible  depends  on  the  form  of 
V.(p),  and  more  particixlar3y on  V'(p,i^).  If  V.  is  independent  of  a,  p,  the 
polar  coordinates  in  the  p- coordinates,  vre  have, 

(30)  0Jp)   =  0^(p)  Yj^Ja,p); 

where  Y.   are  normalized  spherical  harmonics.  Then  if  V  is  independent  of 
(a,p),  j^  =  j^;  if  V  ru  cos  a,  J^  =  Jf  ±  1.  etc. 

(iv)  Overall,  the  effect  of  molecular  rotation  on  the  vibrational  motion  will 
be  the  addition  of  rotational  fine  structure  to  the  energy  level  diagrams . 
However,  this  does  not  invalidate  any  of  the  present  arguments  regarding  the 
character  and  magnitude  of  quantum  jumps,  because  of  the  selection  rules  that 
will  be  enforced  by  a  'reasonable'  form  of  V.  One  will  presumably  have  to 
make  a  number  of  calculations  like  the  present  one,  rather  than  take  combin- 
ations of  possible  transition  matrix  elements  within  each  calculation. 

9-   Discussion 

In  the  preceding  analysis  the  problem  of  calculating  reaction  rates  has 
been  examined;  we  see  that  it  is  a  very  difficult  problem.  Among  other  things, 
we  have  seen  that  the  customary  statistical  theory  is  inadequate  in  many 
respects.  Yet  it  is  perfectly  reasonable  for  a  chemist,  who  is  familiar  with 
the  statistical  theory  and  knows  its  successes,  to  say  that  I  admit  that  I 
cannot  yet  make  qviantitative  calculations  of  reaction  rates,  yet  I  object  to  the 
statistical  theory  which  does  make  quantitative  predictions.   Or,  to  put  it 
another  way,  why  does  the  (supposedly  inadequate)  statistical  theory  give 
sensible  results? 

The  answer  to  this  has  several  aspects.  First,  the  statistical  theory 
works  at  the  phenomenological  level  of  explaining  qualitatively  the  meaning 
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of  the  two  parameters  E  ,  and  P„  in  the  expression  (2)  for  the  reaction  rate 
k_  (T),  and  consequently  of  correlating  E    and  P.  values  for  different  re- 
actions.  It  is  clear  that  once  one  knows  the  E-  and  P- values  for  a  number  of 
reactions,  one  can  go  a  fair  distance  towards  interpolating  rates  of  related 
processes,  or  even  making  moderate  extrapolations;  but  if  there  is  no  information 
that  bears  on  a  given  sort  of  reaction,  or  if  simple  correlations  do  not  give  the 
correct  result,  it  is  necessary  to  go  to  some  deeper  theory,  such  as  the  present 
one. 

Furthermore,  once  we  have  a  potential  energy  surface  for  a  given  reaction, 
it  will  be  possible  to  use  both  the  classical,  statistical  theory  and  the 
present  quantum-mechanical  theory  to  predict  E-  and  P- values.  These  values 
will  evidently  be  different  on  account  of  the  different  models  used,  but  to 
make  a  comparison  between  the  two  theories  that  can  be  compared  with  experiment 
we  should  need  very  much  better  potential  energy  functions  than  are  available 
at  present  (cf.  Section  3)- 

In  conclusion,  one  can  paraphrase  the  conclusion  of  a  related  piece  of  work'--^ 
At  present  the  chief  obstacle  to  the  quantitative  calculation  of  reaction  rates 
is  the  lack  of  detailed  and  accurate  potential  energy  functions.  Once  this  is 
overcome  for  specific  processes,  it  should  be  possible  to  make  rather  accurate 
and  detailed  calculations  of  reaction  rates.  This  will  call  for  considerable 
effort,  though  very  much  less  than  is  involved  in  obtaining  potential  energy 
surfaces.  The  statistical  theory  will  certainly  remain  useful  at  the  empirical, 
level  outlined  above,  though  for  certain  elementary  and  crucial  reactions  such 
as  (l6a,b),  the  present  method  may  well  be  the  only  way  of  getting  reliable 
theoretical  answers. 
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